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1.  Introduction 


This  preliminary  report  is  concerned  with  an  approximate 
description  of  the  two  dimensional  vortical  motion  in  a  thin 
spherical  layer  of  a  gravitating,  incompressible,  inviscid  fluid. 
We  suppose  that  the  fluid  is  contained  between  an  inner  rigid 
ball  and  an  outer  rigid  concentric  sphere  and  that  the  fluid, 
along  with  the  inner  spherical  surface,  has  been  rotating  with 
constant  angular  velocity  about  a  polar  axis  until  a  reference 
time  t  =  0.  Thereafter  we  suppose  that  the  motion  is.  due  to  the 
gravitational  attraction  of  the  rotating  ball  and  to  the  sudden 
creation  at  t  =  0  of  concentrated  vortices  with  axis  normal  to 
the  spheres. 

The  general  objective  is  to  discuss  the  subsequent  motion 
of  these  isolated  vortices,  including  the  case  for  which  cu  =  0. 
However,  in  this  report  we  confine  the  discussion  primarily  to 
the  motion  of  a  single  concentrated  vortex  in  the  northern  hemi¬ 
sphere  when  certain  equatorial  boundary  conditions  are  imposed. 
The  differential  equations  for  the  motion  of  such  a  vortex  appear 
below;  and  they  are  followed  by  an  approximate  description  of  its 
geometric  path  for  a  period  of  time  of  such  duration  that  i"t  is 
not  necessary  to  use  a  computer  for  the  numerical  integration  of 
the  equations. 

There  are  many  physical  phenomena  whi'-'h  are  related  to  the 
motion  of  vortices  in  a  rotating  fluid;  and  the  analysis  of  this 
kind  of  motion  poses  various  inherently  interesting  mathematical 
problems.  Here,  the  investigation  was  undertaken  with  the  hope 


that  the  results  could  be  used  to  approximate  roughly  the  path 
of  a  hurricane  generated  in  the  northern  hemisphere. 


2.  Formulation 

L  t  a  and  a+h  be  the  radii  of  two  rigid  concentric  spheres 
which  contain  a  spherical  layer  of  fluid.  Suppose  that  h  >  0  is 
small  compared  with  a;  and  that  the  inner  spnere  rotates  about 
its  polar  axis  with  constant  angular  velocity  oj.  Let  the  motion 
of  the  fluid  be  referred  to  the  inner  sphere.  Let  p  denote  the 
distance  of  a  fluid  particle  from  the  center  of  the  sphere  while 
‘f’  and  ©  respectively  denote  its  longitude  and  colatitude.  In 
terms  'f  these  coordinates  the  velocity  of  a  particle  relative 
■to  the  inner  rotating  sphere  is  defined  by  the  components 


u  -  (p  sin©)  ^  =  tangential  component  toward  the  east; 
d© 

v  “  "P  cTE-  =  tangential  component  toward  the  north; 


w 


=  radial  component. 


If  the  only  body  force  acting  is  that  due  to  the  gravita¬ 
tional  attraction  of  the  solid  ball  defined  by  the  inner  sphere, 
then  the  basic  hydrodynamical  equations  which  define  the  motion 
of  the  fluid  are  the  continuity  equation 

_1  bp1'  w  +  1 

p  dp  sin  © 

and  the  momentum  equations 


du  d(v  sin  0  ) 
Jl  ~  bt 


2 


du  uv  cot  0 


dt  p 


dV  VW 

HT  4'  “  + 

2 

dw 

"3T  -  p 


+  ™  + 


u^cot  0 


+ 


—  -  2cjdu 

P 


Pcuw  sin  6  - 


2cdu  cos  0 


sin  0  =  -g 


2o)V  cos  0  = 


1  dp  . 

6jP 

1  dp 

"  6o  ' 


1 


5  p  sin  0 

o' 


dp 
d£  ’ 


In  these  equations  the  differential  operator  with  respect  to  the 
time  means 

d_dL  u  d  v  d  d 

13X  "  "5Z  p  sin”0  “  p  +  w  Sp  • 


The  symbol  6  denotes  the  constant  density  of  the  fluid;  p  denotes 
the  pressure;  and  gp  denotes  the  gravitational  potential  modified 
by  the  centrifugal  effect  of  rotation. 

Since  h  is  small  compared  with  the  large  radius  a  of  the 
inner  sDhere;  and  since 


w(  <j>,  0,  a+h,  t )  =  w(|>,0,a,t)  =  0 


let  us  assume  that  the  radial  velocity  and  the  radial  variation 
of  u  and  v  can  be  neglected.  Lt  1  us  also  assume  that  the  motion 
is  such  that  the  nonlinear  terms  in  the  tangential  momentum  equa¬ 
tions  can  be  neglected;  and  that  the  radial  momentum  equation  can 
be  replaced  with  the  hydrostatic  law 

p(<}>,P,p,t)  =  p(;J),e,a,t)  -  g6o(p-a)  . 

Under  these  assumptions  which  characterize  the  linear  shallow 
water  theory;  and  with  the  notation 
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u(<M,t)  =  u{  (f>,  0,  a,  t )  , 
v( <|>, e,  t )  =  v(ij>,e,a,  t)  , 
P(<M*  t)  =  p(4>,  a,  t)  ; 


an  approximation  to  the  motion  is  determined  by  the  equations 

(2.1) 


du  d(v  sin  6) 

af - - -  0  ; 


(2.2) 

(2.3) 


bu  0  - 

-  2<w  cos  6  =  - 


1  bp  . 
6oa  Bin  e  b$  ’ 


■^r  +  2o)u  cos  6  =  — — 
ot  6  a  ao 

o 


Ihe  elimination  of  p  from  (2.2)  and  (2.3)  and  the  use  of 
(2.1)  leads  to  the  equation 


(2.4) 


S  1 
"St"  a  sin  6 


d(u  sin  6 )  ,  dv 

- 5B -  +  ^ 


2o)V  sin  6 

-  —  u  • 


It  should  be  remarked  that  the  quantity 


C  = 


1 

a  sin  b 


S(u  sin  0)  .  bv 

5? - 


is  the  radial  component  of  vorticity  when  radial  variations  are 
neglected.  An  integration  of  (2.4)  yields 


(2.5) 


a  sin  6 


bu  sin  0  ,  bv1 

— sr —  +  aj 


2o)  r 
a  J 


v  sin  0  dx  =  C(<j),6,0) 


0 
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where  C(|),  6,0)  represents  the  vorticity  prevalent  in  the  fluid 
layer  at  time  t  =  0. 

The  equations  (2.1)  and  (2.5)  are  the  basic  ones  for  the 
motion  under  consideration.  They  can  be  studied  by  introducing 
the  function  ^({>,  6,  t )  and  setting 


(2.6) 


(a  sin  0)  = 


~  1  dip 

U  —  —  -57T 

a  dP 


(2.7) 


d0  ~  1  bif/ 

a  IT  "  v  "  a  sin"?  ' 


in  this  way  (2.1)  is  automatically  satisfied.  Then  if  we  sub¬ 
stitute  (2.6)  and  (2.7)  in  (2.5)  we  find  J  must  satisfy 


(2.8) 


1  a 

.n  0  "3^ 


(sin6)||  +  Vr-  J  ^dt  _  a2C(|>,0,O)  . 

sm  6  eij>  n  v 


The  problem  now  is  to  find  a  function  if/  which  satisfies  certain 
prescribed  conditions  along  the  boundary  of  a  domain  D  on  the 
reference  sphere;  and  satisfies  (2.8)  at  each  point  of  D.  Before 
we  consider  this  problem,  however,  let  us  discuss  briefly  the 
simpler  case  in  which  the  spheres  are  fixed  and  the  contained 
spherical  layer  of  fluid  is  motionless  until  t  =  0. 


3.  The  Case  of  no  Rotation 

If  oo  =  0  the  equation  (2.8)  become: 


(5-1) 


^  (slnC)^  +  ^|g  =  a7U,O0) 
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It  is  interesting  to  observe  that  if  we  use  the  transformation, 
(see  Fig.  1), 


z 


r  e 


^  =  (a  tan  ^ ) e ^ ^ 


x  +  iy  , 


then  (3*1)  can  be  reduced  to  the  ordinary  two  dimensional  poten¬ 
tial  equation.  For  example,  suppose  we  wish  to  analvze  the 
motion  of  a  concentrated  vortex  in  the  zonal  layer  defined  by 

ra  <  6  <  £ 

D  :  ] 

{  0  <  <j)  <  2?r 


on  the  reference  sphere  subject  to  the  requirement  that  along 
the  boundary  the  normal  velocity  v  must  be  zero  along  0  =  a  and 
0=0.  If  we  use  the  above  transformation  the  problem  can  be 
reduced  to  solving 


(3.2) 


(x  -  x1)b(y  -  yx) 


for  the  annular  domain 


a  tan  ^  <  r  <  a  tan  ^  =  q,. 


D  : 


0  <  <  2tt 


with  the  boundary  conditions 

Vitiyk)  =  ^(q. .,<!>)  =  o  . 

In  equation  (3 .2)  q  represents,  the  strength  of  the  vortex;  6 
stands  for  the  Dirac  delta  function;  and  x^,  y^  correspond  to 
the  position  of  the  vortex  on  the  sphere. 
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It  should  be  remarked  that  the  annular  domain  D  is  a 
stereographic  projection  of  the  zone  D  as  Fig.  1  indicates. 

The  above  observatiors  suggest  the  possibility  of  de¬ 
veloping  a  general  theory  for  the  motion  of  concentrated  vortices 
on  a  fixed  sphere.  In  the  course  of  doing  this  the  author  found 
that  it  is  basically  analogous  to  the  Helmholtz-Kirchof f  theory 
for  the  planar  motion  of  rectilinear  vortices. 


4.  Fundamental  Solution.  Concentrated  Vortex  on  a  Rotating 
Sphere 

Let  us  return  now  to  the  equation  (2.8)  which  can  be 
written  in  the  form 


(4.1) 


_  ,  _  1  d  ,  „ .  p  \  d}  ,  1 

t  ^ 

C(<tbG>0)  -  2o >  J  dT  • 

0 


=  a 


This  equation  can  be  analyzed  in  several  different  ways.  One  way 
is  to  introduce  an  appropriate  Green's  function  and  then  present 
the  partial  differential  equation  as  an  integral  equation.  The 
Neumann  expansion  of  the  integral  equation  leads  to 

co  /  ,  \n 


n=l 


ni  rn 


It  can  be  shown  that  the  function  in  this  expansion  must  satisfy 
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Oro 


0>o  =  a^(<{>,e,0)  , 


** n  ■  -2  af  Vi  ; 


and  therefore,  by  virtue  of  the  transformation 


.  0  \  ii> 

z  =  (a  tan  7y)e  T  , 


they  can  be  determined  by  solving  a  sequence  of  ordinary  pot*  n- 
tial  problems. 

Another  way  to  solve  (4.1)  is  to  use  transform  theory  based 
on  finite  Fourier  transforms  in  conjunction  with  Legendre  1 ran^- 
forms.  This  way  seems  most  convenient  for  our  immediatt  aim 
which  is  to  find  the  function  $  associated  with  the  creation  at 
t  =  0  of  a  vortex  of  strength  p  concentrated  at  where 

0  <  0-j-  <  TT  . 


In  other  words,  our  first  object  is  to  solve 


(4.2) 


(sin  G)  It  +  — V-  ^  +  2a) 

sin  4  do  sirp  e 


0 


dG 


dT 


M-5  (  4>  -  <t>q)6(G“  0q) 

sin  44 


for 


0  <  0  <  7 r  ,  0  <  (f>  <  2tt  ; 


that  is,  for  the  whole  sphere  excluding  the  south  pole  whpre 
Q( 4»,  e,  t )  may  be  singular. 

Let  us  assume  that  the  behavior  of  G  at  th*  south  p  >le  is 
such  that 
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sin  6  G ( 4>,  6,  t )  =  0  , 


Le->rr  sin  6  If 


=  const, 


Under  these  assumptions  it  is  known  that  G  can  be  expressed  in 


the  form 


(4.3)  G(<M,t)  = 


J— 

n=l  m=l 


A(0, t) 


„  A  (t)cos  mi 

—  (n-m)l  ™  ? 


+BnIn(t)sin  m(J> 


F^(cos  e ) 


where  the  coefficients  are  the  transforms 


c.u 

A(0,  t)  =  J  G(a,0,t)da  ; 


2tt  tr 


AnIn(t)  ~JJ  G(a,  (3,t)cos  ma  P^(cos  f3  )sin  [3  d(3da  ; 


0  0 
2ir  it 


B  m(t)  =  /  /  G(a,  (3,  t  )sin  ma  A  cos  )  sin  P  d(3da  . 

run 

0  0 

The  function  F^(cos  0)  is  the  associated  Legendre  spherical 
harmonic  of  degree  n  and  order  m.  It  can  be  defined  by 


I^(cos  0)  =  < 


(-l)m( sin  0)ni(J  ~~  7T )  P  (cos  0)  ,  m  <  n 


m  >  n 


and  it  satisfies 
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sin 


+ 


~u  w  ^sin  e)  w  ^cos  6 ) 

2  -1 


n(n^-l)  - 


m 


=  0  . 


c  i  n  £) 


P^(cos  0) 


Notice  that  no  previous  condition  on  G  is  affected  if  we  change 
G  by  a  constant. 

The  transform  A(i9,t)  must  be  continuous  at  0=9^  and 
satisfy 


(sin  6 )  If  -  46 (  ~  ei') 


Since  h(6,t)  must  be  bounded  at  6  =  0,  the  solution  of  the  last 
equation,  excepting  an  additive  constant,  is  uniquely  given  by 

e. 


(4.4) 


A(G)  = 


p  £n  tan 


2  ' 


p  In  tan  ^  , 


0  <  , 


°1  «  ° 


This  function  can  be  expanded  in  che  form 


e 

In  cos  — 


(4.5)  A(G)  = 


on  , , 
u  <  — 


n= 


-  P„(cos  e1))Pn(cos  0) 


n 


where  Pn(  )  =  pn(  )• 

The  transforms  A  (t)  and  B  (t)  must  satisfy 

t 

-n(n+l)Anjn(t )  +  2oum  J  B^frJdT  =  p  cos  mcf^-P^cos  )  , 

0 
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-n(ntl  )Brun(  t )  -  2am  J  AnjTi(x)dT  =  a  sin  pJJ(cos  0± ) 
and  these  yield 

dP^(cos  0-,  ) 

Anm^  =  - 

(4.6) 


0 


Anm(t)  ”  -  — "n(n+ir~  C0E  ”  [+1  ' 


2o)t 

n(ntl)J  ’ 


hpJJ(cos  e1) 

Bnm(t)  = - n(n+l )  sin  ” 


d)  -  2a)t 
.  1  n(  n+1 ) 


The  substitution  of  these  transforms  in  (4.3)  gives 


(4.7)  G(<M,t)  = 


n 


A(e) 

00  (2n+l)  J±-  (n-m)I 


^  '  ^1 +  n['n+l )]  • 


.  cos  m 


Next,  the  addition  theorem 


(^•8)  Pn[cos0cosei  +  sin  6  sinG^  cos  A] 

Pn(cos  0)Pn(cos  e1) 

+  2  ^  Tn+mii  ^(cos  VOcos  0  )cos  NA 

shows  that  we  can  write 


12 


(4.9)  0(4),  e,  t )  = 


n 


U  (2n+l)  } 

"  W  n(n+i)  \ 


A(0) 

cos  0  cos  0. 


+sin  0  sin  0;L  cos(  4  -  ^  +  n 7 \ j 


-Pn(cos  01)Pn(cos  0) 


Hence,  from  (4.5),  we  see  that  G(4,0,t)  can  be  exhibited  in  the 
form 

0. 


(4.10)  g(4>  0, t )  = 


H  ^  (2n-fl) 

4tt  “y  n  ( n+ 1 ) 


n 


-  in  cos 


cos  0  cos  0. 


2u>t  ^ 


+sin  0  sin  01  cos  ( 4  -  4p  + 
-(-l)nPn(cos  0) 


Finally,  by  using 


(4.11) 


we  find  that 


QD 


x -  ( 2n+l )  /  i )np  i v \ 


( 1+x ) 

1  +  £n  2 


(4.12)  G(4,0,t)  = 


H  fl 

"  3r|j? 


|i  t ( ^+1  )  P 

W  n(n+1-)  n 


i  0 

£n  cos  -g-+  £n  cos 


cos  0  cos  0. 


+sin  0  sin  01  cos (4  -  4p  +  n(nVl'j) 
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is  that  when 


1  i 

The  reason  for  retaining  the  constant  £+  |n  cos 
we  do 


G (<M,  t)  =  G((j),e,t;4)1,01) 
is  symmetric  in  6  and  6^. 

Except  for  ( and  the  south  pole,  the  velocity  compo¬ 
nents  are  given  by 

~  1  dG  .  v  =  1  dG 

a  ’  a  sin  6  ^  * 


The  velocity  field  in  the  neighborhood  of  the  concentrated  vortex 
can  be  estimated  by  studying  the  first  term  of  the  expansion  of 
G ( 4>,  0,  t )  in  powers  of  cot,  namely 


(4.13)  G  U,0,t)  = 


a  f  ft' 

-  1  +  2  |n  cos  -£-+  2  |n  cos  ^ 


QD 


4  (2n+1)  p 

tijr  n(  n+1  j  n 


cos  6  cos  6- 


+sin  0  sin  0^  cos(<j)  -  <J>^) 


a 

4tF 


p  ^  1  OP 

•1  -  in  cos  -  in  cos  ^ 


+1  +  |n 


1-  cos  0  cos  0^  -  sin©  sin0^cos((J>-  ^  ) 

5 - 


=  jp  In 

4tt 


2  © 


2ei 


j^tan  -£+tan  -  2  tan  ^  tan  cos  ( <J;  -  ^  ) 


This  gives  the  dominant  part  of  the  velocity  of  a  point  near 
(<l)l>®l)*  It  corresponds  to  the  case  of  no  rotation;  and  if  we 
transform  to  the  z-plane  introduced  above  via 


l4 


z 


we  have 

(4.l4) 


6  i(J> 

=  a  tan  2  e 


Go  "  W  £n  "  Z1^Z  '  Zl^  "  4tF  1 


H  ...  _2 


n  a 


=  Re  «n  (z  -  z,  )  -  tc:  in  a 


c.7T 


TT 


where  zx  is  the  image  of  (<1^,^)  and  Re  P(z)  denotes  the  real 
part  01  p(z).  The  velocity  components  determined  by  C.Q  can  be 
calculated  from 


a  dO 


n  _ 

p(  a  +  zz  ) 


Re 


r\j 


1 _  *Go  =  n(aP+  zz)  Re 

TnT  dJT  2  H1 

4ira  /zz 


These  show  that  as  z  — *■  z-^  the  velocity  approaches  ini  ini ty,  but 
ti'.o  component  In  the  direction  from  z  to  ■i1  is  zero.  Therefore  m 
the  neighborhood  of  the  fluid  rotates  about  the  vortex 

point  with  velocity  of  the  order  of  p  Fro"'  this,  it  can  be 

deduced  that  the  velocity  field  defined  by  G  in  the  neighborhood 
of  behaves  qualitatively  in  the  same  way  as  the  field 

defined  by  G  .  In  other  words,  the  concentrated  vortex  at  {^,6^ 
has  no  tendency  to  move  —  it  is  not  sell  propelling. 

Notice  that  if  the  vortex  is  at  the  north  pole,  i.e., 


6^  =  0,  then 

G((J>,0,t;i|)^,0) 


-  —  \: 

?rrl: 


1 


+ 


in  cos 


P  (cos  0 ) 
11 v 
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G(M»tj  0)  =  Go(4»,  0,  t;  4>1,  0) 


£n  tan 


0 

2 


This  function  is  such  that 


J0  ->0 


sin0  G0  ( 4>,  ©,  t;^,  0)  =  — 


5 •  Concentrated  Vortex  in  Northern  Hemisphere.  Either  Velocity 
Component  Zero  at  the  Equator 

We  proceed  to  find  the  stream  function  for  a  concentrated 
vortex  traveling  in  the  northern  hemisphere  and  subject  to  the 
condition  that  the  normal  velocity  of  the  fluid  at  the  equator  is 
zero.  The  boundary  condition  then  is 


( <J>*  ^  *  t )  - 


a  sm 


or 


,t)  =  0 


5  > t) = 


However,  we  know  that  nothing  is  lost  if  we  impose 

$(  <ib  \  =  0  . 


The  stream  function  for  this  case  can  be  formed  from 
G(  (J>,  0,  t;  if>^,  0^ )  if  we  assume  that  the  principle  of  reflection 
across  the  equator  will  lead  toward  the  desired  result.  In  fact, 
it  is  easy  to  verify  that 

(5*1)  tP  =  G(<}>,0,  t;<J)^,01)  -  G((|),0,t;|i^,Tr-  0-^ ) 
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is  the  desired  function.  We  refrain  from  a  discussion  of  the 
velocity  field  defined  by  (5.1)  in  favor  of  the  next  case  in 
which  the  eastward  velocity 


is  prescribed  at  the  equator. 

We  turn  to  the  problem  of  finding  such  that  with  the 
boundary  condition 


the  equation 
(5.2) 


a  •  ( <tb  TT  >  t )  =  c  =  const. 


1  d  /  „ ^  \  ,  1  d" 


t 


sITTe  w  vGin  e)  w  + 


bip 


- ^4+  2 0)  /  ^dx 

sin^e  ^tj)  d<j> 

m.6(4»  -  4>1)b(©  -  ©1) 

sin  07 


is  satisfied  for 


If  we  assume  that 


(5.3) 


0  _<  t{>  <  2xr  ,  0  <  Q  K  1L  . 

—  2 


4, 


’-•0+ o  Ein  e*n  =  - 


6  2tt 


ai  integration  of  (5«2)  over  the  northern  hemisphere  yield: 

2tt 


(5-^) 


J  ^  ^  t  )d^>  =  q  +  (j,Q  . 


0 


This  shows  tha.  the  prescription  of  TpQ(i j),  J,t)  must  satisfy  the 
compatabi lity  condition  (5*^)»  If  we  wish  to  satisfy 
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V0(b,  \>  t)  =  ac 


then  we  must  have 


b0  =  2rrac  -  g 


As  we  have  seen,  the  condition  (5*3)  is  satisfied  by 


—  G(  4>,  0,  t;  0)  -  -^2.  fn  tan  — 


Consider  the  function 


(5.5)  = 


in  tan  ^ 


+G(  4»,  o,  t;  t^1,  ei)  +  G(  4>,  e,  t;  -  0-l) 


A  short  computation  chows  that 


Hence  (5.5)  satisfies  the  boundary  condition  that  the  eastward 
velocity  at  the  equator  is  constant  and  equal  to  c  if 


(i  =  2mac  -  g  . 

We  have  now  shown  that 


(5.6)  tf'('tb^t)  = 


( 2-rrac  -  g  ) 
2tt 


in  tan 


6 

7 


+G(  {|>,  0,  t;  <|>1,  ©1 )  +  G(<J>,0,  t;  ^,-TT  -  0-j  ) 

is  the  desired  stream  function. 

Let  us  turn  our  attention  to  a  study  of  the  path  of  the 
vortex  at  [  t)>1(  t ),  01(  t )] .  As  in  the  basic  Helmholtz  theory  for 
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the  motion  of  interacting  concentrated  vortices,  the  velocity  of 
the  vortex  at  [ $  ( t ), 0  ( t ) ]  is  equal  to  the  field  velocity  deter¬ 
mined  at  [(^,6-jJ  by  all  of  the  other  vortices.  That  is,  the 
equations  for  the  motion  of  the  vortex  are 


^5 

a  sin  0. 


3 


i2™LZ*l  {n  tan  § 

[_+a ( 4>, e,  -  ep 


and 


1=^ 

6=6. 


a 

a  sin  0 


—  (  |{),  0  ,  t ^ ,  7T  -  )J 


0=0- 


In  expanded  form  these  are 


*  1  ^ 1  2rrac  -  u 

(5*7)  ^  =  -fir-  = - £ - z —  + 

1  at  2rra  sin^01 


_ U 

2  2  ^l 

Bma  cos 


u  cos  0X  ^00.  (2n+1) 

”  2  ‘ — T  n  ( n+ 1 ) 

4ma  n=l 


n 


-cos'1 0. 


+sin‘*0n  cos 


2cut 


n (n+TJ 


and 

d0, 

(5»8)  ^  “  -tf- 


•  1  +  cos 


2a)t 

n(  n+1 )_ 


u  sin  tp  CD  (gn+1) 
4ira'r‘  fcr  nTii+rj 


n 


-cos'-©. 


+sin‘-01  cos  n^-fT j 


•  sin 


2cot 
n  ( n+ 1 ) 
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It  seems  that  numerical  integration  is  necessary  in  order 
to  trace  the  path  of  the  vortex  for  moderately  large  or  large 
values  of  oot.  Our  intention  is  to  pursue  this  in  a  later  report; 
however,  as  a  preliminary  step,  we  can  assume  that  oot  is  so  small 
that  a  sufficiently  good  approximation  to  the  path  can  be  found 

e 

by  expanding  ^  and  in  powers  of  oot  and  retaining  those  of 
degree  no  higher  than  the  second.  In  addition  let  us  suppose 
hereafter  that  c  =  0. 

If  each  term  in  (5*7)  and  (5*8)  is  expanded  in  powers  of 
oot;  and  if  we  retain  powers  no  higher  than  the  second  we  find 


(5.9)  = - — n —  +  - ± - rr- 

2r ra  sin  0,  D  2  2  P1 

1  Hma  cos  -7T- 

d 


llCOS0lr®  ( 2n+l )  1 

— — 4-t-  nTn+TT  < 


(a)t ) 

,  n^(n+l)‘ 


P^l-cos  2G1) 


P^-cos  2 0±) 


+2  sin^'G-^.  P”  ( -  cos  2 0^ ) 


.  d(ojt)sin  0  oo  /p  ,  . 

(5‘10)  =  ~  S^Pn(~e0Sa9^  ’ 


The  series  which  appear  in  the  last  equations  can  be  summed  and 
they  lead  to 
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(5.12) 


jioot  in  cos  0^ 
2jra~  jl  -  cos261 


In  order  to  express  ^  and  6^  in  powers  of  t  explicitly  we  need 
to  expand  cos  0^  in  powers  of  t  and  insert  the  expansion  in 
(5.11)  and  (5.12).  If  we  do  this  and  retain  only  pov/ers  of  t 
less  than  the  third  we  find 


(5.13)  - - £ 


u 


8  sin 


2^ 


Tra  cos  0 


o 


+  (o>t )‘ 


I 


cos  0 


0  }?  in  \ 


d* 


ll ( 1  -  2  cos  0  )  in  cos  0 
^ v  o  o 

- - - pr - 

8ira  cd  sin  cos  0q 


J 


l-uut  in  cos  0 

(5.14)  0  =  - ^ - - 

2rra  sin  0Q 

where  0  is  the  colatitude  of  the  initial  position  of  the  vortex, 
o 

The  approximate  equations  (5*15)  and  (5»l4)  show  the 
following  if  \i,  the  strength  of  the  vortex,  is  positive;  and  0q 
is  greater  than  v/j>.  When  cut  is  very  small  the  vortex  moves 
along  a  path  which  is  almost  parabolic.  Since  -a.6^  is  the  vel¬ 
ocity  to  the  north,  equation  (5.14)  shows  that  the  rate  of  change 
of  latitude  is  zero  at  t  =  0  and  the  latitude  increases  as  t 
increases.  Equation  (5.15)  shows  that  the  vortex  starts  with  a 
westward  motion  but  after  a  while  the  positive  term 
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n(oot)' 


wac  cos  6 


cos  6 


0  J  in  X  ,  'i{1-  2  C0E  V*11  co°  °o 

7"  d?s  +  - rr 


0 


1  -  A 


2  4  o 

Brra  a)  sin  -p-  cos  0 


becomes  dominant  and  the  vortex  then  moves  to  the  east.  The 
following  figure  shows  a  qualitative  sketch  of  the  path. 


Figure  2 


This  sketch  exhibits  some  of  the  characteristics  cf  the  observed 
paths  of  hurricanes  generated  in  the  Caribbean  sea  or  some  other 
southern  part  of  the  North  Atlantic  ocean. 
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